Point- Vortex Simulations Reveal Universality Class in Growth of 2D Turbulent 

Mixing Layers 
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A central but controversial issue in free turbulent shear flows has been the universality (or other- 
wise) of their growth rates. We resolve this issue here in the special case of a temporal 2D mixing 
layer in a point vortex gas by extensive high-precision numerical simulations, utilizing for the first 
time a powerful ensemble-averaging strategy. The simulations show that the momentum thickness 
of such a mixing layer grows at the universal asymptotic rate of 0.0167(±0. 00017) times the velocity 
differential across the layer over a wide range of initial conditions, often after very long transients. 



The quest for universality in turbulent flow goes back 
to Reynolds [l[ , and includes the discovery of the log law 
in wall-bounded flows by Prandtl and Karman Q and the 
£-5/3 S p ec t rum by Kolmogorov [3] . In turbulent shear 
flows, a unique 'equilibrium' state independent of the de- 
tailed initial conditions has often been postulated [|| , but 
remains a controversial issue in wall-bounded 0,0] as well 
as in free shear flows (e.g. wakes @, H], mixing layers 0- 
[Hj]). The asymptotic nature of such postulated univer- 
salities, valid only at high Reynolds numbers (Re — > oo) 
and / or far downstream (x — > oo), makes it difficult to 
be certain that the final state has been reached in flows 
with extremely long relaxation times 0, [13} . 

Mixing layers have been widely studied by experiments 
[|,[M[ll r LVortex simulations jlMl , DNS [jlOl^ 
and LES [IJ, [23[ . The dimensionless growth rate based 
on momentum thickness varies from 0.014 to 0.022 across 
experiments [12j | and a similar scatter is observed among 
simulations as well. This has led to suggestions that there 
may be no universal growth rate independent of initial 
conditions [IS Gil- 

We report here results of extensive high-precision sim- 
ulations of a 2D temporal turbulent mixing layer, for a 
wide class of initial conditions involving random and pe- 
riodic normal displacements of a linear, equally spaced 
row of point vortices (separation distance I, Figure [J) 
at the initial instant (t = 0). This can be seen as a 
direct 'molecular dynamics' solution of a class of initial 
value problems we may pose in the spirit of the statisti- 
cal mechanics of a point-vortex gas, formulated first by 
Onsager [24| (see [25[ for a recent review) . Although the 
flow considered is thus 2D instantaneously (and hence 
also in the mean), it is not irrelevant to a plane Navier- 
Stokes mixing layer (3D instantaneously but 2D in the 
mean). This is because experiments indicate that a 2D 
turbulence field provides a reasonable representation be- 
fore the occurrence of the mixing transition [j| I2"6rl28j ] . 
and even later the coherent structures in the flow remain 
quasi- two-dimensional [29| . 

Point vortex simulations of a temporal mixing layer 



Periodic boundary 
conditions 




FIG. 1. Present point vortex model of the canonical tempo- 
ral mixing layer, showing initial perturbation in the form of 
normal displacement of the vortices at t — with amplitude 
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were pioneered by Delcourt and Brown [l5[ and Aref and 
Siggia |l6[, both of whom used cloud-in-cell algorithms, 
and by Acton [l7| ■ We revisit the problem using a differ- 
ent algorithm and the much more powerful computing re- 
sources now available, as the issue about universality still 
remains central to the subject. The present simulations 
use double-precision 4th order Runge-Kutta integration 
to track each individual vortex, and (for the first time) 
provide averages over ensembles, involving upto 108 real- 
izations. These substantial improvements in the numer- 
ics proved crucial for the conclusions we draw here. 

We consider the limiting solution as N ~ > oo, where 
N = L/l is the number of vortices in the streamwise 
domain L (Figure [lj. With periodic boundary conditions 
imposed at x = 0, L, the equations governing the motion 
of the vortices [17| are 
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where T is the (common) strength of each vortex and 
(xi,yi) is the location of vortex i. 

In the first instance the initial displacement of vortex 
i in the array is drawn from a prescribed distribution of 
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random numbers with zero mean and a measure of disper- 
sion denoted by a ; their subsequent motion is obtained 
by integrating (fT] , [5]). The solution for 5 , a measure of 
the thickness of the mixing layer at time t, takes the form 

where F\ is some (unprescribed) function. Several ex- 
ploratory simulations revealed an initial 'transient' solu- 
tion of (P , [5]) that is independent of N at large N . If 
further this solution evolves to a state independent of the 
initial conditions for sufficiently large tAU/l , the depen- 
dence on a/! in (3) will disappear, yielding 

5/L = F 2 (tAU/l), tAU/l » 1 (4) 

where F 2 is an appropriate limit of F± . In what may be 
called the 'outer' limit, tAU/L — O(l), the solution may 
be expected to be of the form 

S/L = F 3 {tAU/L) (5) 

If we now postulate an overlap between solutions Q and 
(J5]) the only possibility is that 

6 = AtAU + B (6) 

where A and B are independent of t. Equation ^ de- 
scribes an 'equilibrium' range in the problem. A univer- 
sality class is defined here as a set of initial conditions 
a/l for which A is independent of a/l in the simultane- 
ous limits tAU/l —> oo and tAU/L — > 0, in the spirit of 
matched asymptotic expansions |3fJf . 

The concept of (fluid-dynamical) equilibrium has been 
defined variously [H, Hl|. We define equilibrium as a 
state in which the mean velocity field and the Reynolds 
shear stress both exhibit self-similarity with the same 
(time-dependent) length and velocity (alternatively time) 
scales. (The Reynolds equation of momentum, see (O 
below, then implies ([S]).) 

In the present code, with the adopted time step of 0.1 
in l/AU, the distance a vortex moves during any time 
step does not exceed that to its nearest neighbor, and 
is almost always at least an order of magnitude less. A 
reduction in time step by a factor of 4 did not materially 
affect the results. After every 100 time steps the x and y 
components of velocity (u, v) are computed on a grid of 
20L points in x and 200 points in y using the Biot-Savart 
relation, and x-averaged quantities like the mean- velocity 
U(y,t) = (l/L) dx[u(x,y,t)] are computed. The so- 
called momentum thickness, defined as [32j j 
. e(t) = JZ o dy[0.25-(U(y,t)/AUf] 
is often preferred here as a measure of layer thickness 
as the integral makes it more robust than the vorticity 
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FIG. 2. A. Effect of ensemble averaging. B. Variation of RMS 
of relative departure of 8 from (8) with n (N — 3200) and N 
(n = 32). C. Growth rates of two different ensembles(n = 10) 
with same class of initial conditions and domain size; small 
ensembles of 10 members and fits over narrow data ranges 
(ti/t\ = 2, where t\ and ti are beginning and end of linear 
regime used to assess the growth rate) can result in large vari- 
ations (14% in this case) in equilibrium growth rates. Note 
that fe/ii and variation seen here are comparable to those 
presented by Balaras et alfiH 

thickness and other point-based measures. However at 
tAU/l < 10, can be misleading because of large over- 
shoots in the mean velocity profile, and S (defined as the 
separation in y between the vortices at extreme y posi- 
tions) would be a more appropriate choice. (This explains 
its use in Regime I below.) 

The accuracy of the algorithm has been assessed in 
several ways. The Hamiltonian, given for a cloud of N 
point vortices by 

" H = -(r 2 /2^)E 2<J ln(| r,-rj | /L),r, = (x uyi ), 
exhibits a maximum deviation (at tAU/l — 2500) of 9 x 
10~ 6 of its initial value for N — 3200. The first moments 
of the vorticity distribution about the x- and y-axes are 
conserved to within 10~ 16 and 3 x 10~ 13 times I, and 
the second moment to within 1.3 x 10~ 9 of its initial 
value. Another check is provided by mean momentum 
balance. For the present model the Reynolds-averaged 
Navier Stokes equation for ^-momentum simplifies to 

where we have used the Reynolds decomposition u — 
U + v! etc., and uj = uJ + lj' is the total vorticity (note 
that V = 0). We find that the integral of 0, in what 
is defined as Regime II below, is satisfied to better than 
0.05 % (in simulations using 1600 vortices). These num- 
bers show that the current computations are substan- 
tially more accurate than any previous work. 

Finally, we note that averaging over a sufficiently large 
domain (analogous to long time averages in spatial mix- 
ing layer experiments) is essential to obtain accurate es- 
timates of growth rate in temporal simulations. Averag- 
ing over a large ensemble of different realizations is an 
equivalent alternative (Figure UK)- In the current sim- 
ulations, the RMS fluctuation from the mean decreases 
with ensemble size n like n~? (Figure 03 )> whereas the 
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FIG. 3. Composite diagram showing effect of initial condi- 
tions and domain size on the evolution of the mixing layer. 
Note use of 5 and 8 in different parts of the diagram, and 
change in the abscissa from tAU/l with a logarithmic scale 
upto 500, linear scale between 500 and 1500, and a switch to 
tAU/L thereafter. Appropriate changes have been made on 
both abscissa and ordinate to ensure that the evolution curve 
should go smoothly from one regime to the next. Inset on top 
left gives variation of initial transient with the amplitude of 
the initial vortex displacement. Two insets on the right give 
pictures of the configuration of the vortices at tAU/L — 1 
(upper) and 4 (lower). 

computational effort increases like n. On the other hand 
the statistical fluctuations again decrease with number 
of vortices as , whereas the computational effort in- 
creases more rapidly like N 2 . Once N is sufficiently large 
the former approach is thus computationally more eco- 
nomical, and is adopted here. 

Figure [3] illustrates the evolution of the flow through 
results from selected simulations, with different domain 
sizes and initial conditions (a/l ranging from 10~ 4 to 10). 
It is seen that the variation of 6 exhibits three distinct 
temporal regimes. 

In Regime I, different initial conditions lead to widely 
different growth histories of S from t = to a value of 
tAU/l that depends on a/l. The duration ti of the tran- 
sient (defined as the time at which the departure of 6 or 
8 from the respective line of best fit in Regime II drops 
below 10%) varies by an order of magnitude, and its de- 
pendence on a/l exhibits a minimum. On the lower side 
of this minimum the disturbance is too small to quickly 
trigger turbulence, and on the other side it is too large 
to die down quickly. There are therefore optimal 'trips' 
that lead to shortest transients. 

In Regime II, which is between Regime I and tAU/ L ~ 
1, growth is linear and independent of the initial con- 
ditions as well as domain size (this claim will be more 
elaborately supported below); in other words, the mix- 
ing layer is now in 'equilibrium'. 

In Regime III, beginning at tAU/L ~ 1, the effects 
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FIG. 4. Growth of momentum thickness for different simula- 
tions, shown with t/-axis displaced. The universal equilibrium 
rate ([Sj) is shown in dashed lines, and is obtained by best fit 
to the average of 108 simulations with N = 3200, a/l = 0.05 
(Rl).Note change in scales on both axes, each by the same 
factor, for tAU /I > 1000. For tick marks see footnote on 
Table 1. Rl is shown with two thin lines within ±1% of ([8]) 

of finite domain size become noticeable and the scaling 
length changes over to L. This is due to the small num- 
ber of coherent structures governing the dynamics in this 
regime (see insets in Figure[3]). At lower tAU/L the effect 
is characterized by larger statistical fluctuations. These 
can in principle be reduced by averaging over a larger 
ensemble, but at tAU/L ~ 4 there is only one structure 
left in the domain, with no further opportunity to amal- 
gamate with others and grow. Instead the structure just 
rotates about its own axis , resulting in oscillatory 9 with 
a stationary mean. A rough spatial analog of this effect 
of finite domain size is found in flow experiments carried 
out in smaller wind tunnels [331 ] . 

We now return to more detailed results on equilibrium 
linear growth in Regime II. In order to test whether 
the growth rate is universal, a large set of simulations 
with widely different initial conditions and domain sizes 
have been performed (with each case averaged over a 
large ensemble : details are in Table HJ . The initial vor- 
tex displacement distributions include the bi-modal type, 
in the form of sums of symmetric and asymmetric dis- 
placed Gaussians (respectively BM1, BM2). The respec- 
tive growth histories in Regime II are shown in Figure 
|U We take as reference, the best fit value for the largest 
ensemble simulated here, Rl (n — 108), 



— = 0.0167AJ7 (8) 
at 

which extends over more than a decade in tAU/l. RMS 
deviations from ([5]) are listed in Table HI 

Periodic excitation needs a separate, detailed discus- 
sion, so we include here only one case to show that, 
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TABLE I. 



Code TV 



Initial 
conditions 



Best fit® 
dS ~ 

d(tAU) 



RMS of^ 

relative deviation 
from ^ 



Rl 3200 108 uniform random 0.01667 0.00676 
g/l = 5x 10~ 2 



R2 10000 11 



uniform random 0.01662 
a/l = lx 10 _1 



0.02260 



R3 1600 64 



uniform random 0.01666 
a/l = 10" 8 



0.00835 



Gl 1600 64 



Gaussian^ 
a/l = 5~ 



0.01676 



0.01845 



BM1 1600 



64 Bi-modaF 
ai/l = 1 x io" 1 
a 2 /l = 1 X 10" 1 
d/l = 6x 10 _1 



0.01654 



0.01282 



BM2 1600 



~64 BTmodaF 
ai/l = 1 x io -2 
02/7 = 2 x 10" 2 
d/l = 4x 10" 2 



0.01665 



0.01367 



PI 3200 32 Periodic forcing 0.01667 0.00697 
a/l = Ax 10 _1 
X/l = 100 

a Based on data beginning from the tick mark on each simulation 

in Figure 4 till the end of respective simulation. 
b Defined as minimum with respect to B of 



-E 



, in - no. of data points 



-(0.0167tA!7+-B) ' 
0.0167tAl/+S 
Standard deviation. 
cri, erg - Standard deviations of two Gaussians separated by d. 
0.1% uniform random noise added to generate different 
realizations; A - Wavelength of periodic forcing 



after mimicking experimentally observed behaviour [lOj 
when wavelength is chosen as the relevant length scale, 
the layer eventually grows at the same universal rate (|8|) . 
The time taken to attain universality corresponds to over 
twice the length of the experimental test section in the 



30 Hz case in Figure 13 of Ref. [lOj. It can be concluded 
from Figure |4] that there is a regime of linear growth, 
with a growth rate that is universal to within 1%, for the 
wide class of initial conditions considered here. 

The present simulations yield growth rates of the same 
order as in experiments. This implies that the dominant 
mechanism in the growth of the momentum thickness in 
the mixing layer must be just the kinematics of the Biot- 
Savart relation, which includes the emergence of chaos 
even in a few- vortex system [34T |. Scepticism about uni- 
versality in the growth of (real) mixing layers [lol | stems 
from their known sensitivity to various factors, but an 
important additional factor suggested by the present sim- 
ulations is the inadequacy of flow-development length in 
laboratory set-ups. Such a length is the spatial analogue 
of the long transient noted in some of the present sim- 
ulations. Similarly, as seen in Figure [2jU, short compu- 
tational domains, fits made over small values of £2/^1 
and small ensembles account for the observed variation in 
growth rates. Evidence for universality in the present sim- 
ulations would be weak unless L/l is in excess of 10 3 , av- 
erages are struck over sufficiently large ensembles and the 
equilibrium regime is sufficiently long (Figure [2]). Simi- 
lar factors can account for claims against universality in 



other types of simulation [11]. It is anyway of consid- 



erable interest that there exists at least one prototypical 
turbulent shear flow with three distinct regimes including 
one corresponding to a universal equilibrium state. 
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